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Abstract. Using ideas from algebraic topology and statistical 
mechanics, we generalize Kirchhoff 's network and matrix-tree the- 
orems to finite CW complexes of arbitrary dimension. As an ap- 
plication, we give a formula expressing Reidcmeistcr torsion as an 
enumeration of higher dimensional spanning trees. 



Gustav Kirchhoff's results on electrical networks, which predate 
Maxwell's theory of electromagnet ism, are a product of the mid-19 th 
century [Kil], [Ki2]. Kirchhoff 's network theorem states that in any 
resistive network there is a unique current satisfying Ohm's law and 
Kirchhoff's current and voltage laws, and furthermore this current can 
be explicitly computed. The first complete treatment of the network 
theorem is attributed to Hermann Weyl [W] in 1923. By the mid-20 th 
century, algebraic topology provided key ideas leading to a simple and 
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elegant proof [E],[R],[NS]. A companion result is Kirchhoff's matrix- 
tree theorem which gives a formula for the number of spanning trees 
in a finite connected graph (see [Mo] for a history of this result). This 
paper is an outgrowth of our investigations on the interplay between 
algebraic topology and statistical mechanics [CKS1], [CKS2], [CKS3]. 
Our aim is to generalize Kirchhoff's results to higher dimensions, as 
well as to connect these results to the theory of Reidemeister torsion. 

A high dimensional network theorem. Suppose X is a finite con- 
nected CW complex of dimension d. Let Cj(X;R) denote the cellular 
chain complex of X with real coefficients and the standard inner prod- 
uct ( , ) for which the set of j-cells, denoted Xj, is an ortho normal 
basis. In what follows we fix a function r: X& — > R + ; the value of r 
at a d-cell b is considered to be the resistance of b. Define a linear 
transformation R: Cd(X;R) —> Cd(X;R) by mapping a <i-cell b to r b b 
and extending linearly. Let Bd-i(X;R) C Cd-i(X;R) be the vector 
subspace of (d — l)-boundaries and let Zd(X; R) C Cd(X; R) be the 
vector subspace of c?-cycles. 

Definition 1.1. A network problem for X consists of a choice of p G 
Bd-i(X;M) and q G Zd(X;R), respectively called (d— 1) -boundary cur- 
rent and d-cycle voltage. 1 A solution consists of V, J G Cd{X; R) such 
that 



(1) V = RJ , (Ohm's law) 

(2) dJ = p , (current law) 

(3) (V, z) = (q, z) , for all z G Z d (X) . (voltage law) 

To see why a solution exists, define a modified inner product ( , ) R 
on C d (X; R) by 

(b,b') R =(Rb,b') 

for b, b' G X d . Let 

d R :C d ^(X;R)->C d (X;R) 

denote the formal adjoint to d using the standard inner product on 
C d -i(X; M) and the modified inner product on C d (X; K). Let B R (X; R) 
be the image of d R and note that B R (X; R) is the orthogonal comple- 
ment to Zd{X; R) in Cd{X; R) with respect to the modified inner prod- 
uct. Elementary linear algebra implies d: B R (X;R) — > Bd-i(X;R) 

^^Whcn d = 1, in the terminology of Roth [R], p is a node current and q is a mesh 
voltage, each arising from an external source. Bollobas [B. p. 41] only considers the 
case when q = and p is of the form pii +pjj for a pair of distinct vertices i and j. 
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is an isomorphism. Consequently, there is a unique Jo G B\{X; R) 
such that <9J = p. Set V = R3 . Then (V ,z) = (Jq,z)r = for 
all z G Z d (X;M.). Let Ji be the orthogonal projection of R~ l q onto 
Z<i(X; R) in the modified inner product, and set x = Ji — R~ x q. Then 
(i2Ji - g, 2) = (x, = for all z G Z d (X; R). Set Vi = R3 V Then 
J := Jo + Ji and V := Vo + Vi solve the network problem. It is 
straightforward to show that this solution is unique. 

The above solution to the network problem uses the orthogonal pro- 
jection of Cd{X; R) onto Z d (X] R) in the modified inner product. In 
the classical case d = 1, Kirchhoff gave a formula expressing the or- 
thogonal projection as a weighted sum indexed over the set of spanning 
trees of X. To get an explicit formula in higher dimensions we will need 
a notion of spanning tree. 

Definition 1.2. Assume as above that X is a connected finite CW 
complex of dimension d. A spanning tree for X is a subcomplex T such 
that 

• H d (T;Z) = 0, 

• /3d-i(T) = j3d-i(X), where 0k(X) denotes the k-th Betti num- 
ber, 

• X^" 1 ) C T, where X^ is the fc-skeletion of X. 

Remark 1.3. We will show in the next section such spanning trees exist. 
The reader should be aware that the literature contains sundry notions 
of "high dimensional spanning tree." 2 Note that when d — 1, our 
definition reduces to the classical notion of spanning tree. 

Definition 1.4. For a spanning tree T of X, define a linear transfor- 
mation 

f:C d (X;R)^Z d (X;R) 

as follows: Let b be a <i-cell. If b is contained in T then we set T(b) = 0. 
Otherwise, note that H d (T U 6; Z) = Z d (T U b; Z) is infinite cyclic. 
Let c be a generator. Set t b = (c,b) (this is always non-zero). Then 
T(b) := c/U, is a real d-cycle of X. It is easy to see that T(b) is 
independent of the choice of c. 



2 For example, [DKM] replaces condition (2) with the requirement that the re- 
duced Betti number j3d-\{T) is trivial. This implies Pd-\{X) is trivial as well, so 
the notion of spanning tree in [DKM] does not apply to a general finite complex 
X. Sec [P] for a detailed discussion of the various notions. 
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Let 6t denote the order of the torsion subgroup of Ha-i(T] Z) and 
define the weight of T to be the positive real number 

«T := 01 n r^ 1 . 

b£T d 

Theorem A (Higher Projection Formula). With respect to the modified 
inner product ( , the orthogonal projection Cd(X;R) —> Zd(X;R) 
is given by 

T 

where the sum is over all spanning trees, and A = Y1t w t- 

Let d* : C d -i(X; R) — >■ C d (X; R) be the formal adjoint to the bound- 
ary operator with respect to the standard inner product. Define B d (X, R) 
to be the image of d*. Then we have 

Addendum B (Higher Network Theorem). Given a vector V G Cd{X; R), 
there is only one vector z G Zd(X\R) such that V — Rz G B d (X, R). 
Furthermore, for each d-cell b, we have 

T 

Remark 1.5. In classical network terminology {d = 1), (V, b) is the 
voltage source on branch b and (z, b) is the current resulting in branch 
b (see [R],[NS]). 

A high dimensional matrix-tree theorem. The classical matrix- 
tree theorem enumerates the number of spanning trees of a graph. In 
higher dimensions, the best we can achieve is an expression for Q\. 
Observe that Bd-i(X; R) is an invariant subspace of the operator 

dd* R : C d ^(X;R)^C d ^(X;R). 

Let 

C R : B^X-R) ^ B^X-R) 
denote the associated restriction. 

Theorem C (Higher Weighted Matrix- Tree Theorem). We have 

det C R = 7x w T , 

T 

where the sum is indexed over all spanning trees of X , and the normal- 
izing factor is given by 

fix 
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where fix 6 N is the square of the covolume of the lattice Bd-i(X; Z) C 
Bd-i(X; R) with respect to the restriction of the standard inner product 
ofCd-i(X; R) andOx is the order of the torsion subgroup of H<i-i(X; Z). 

The unweighted case when r : X^ — > R + is constant with value 1 is 
worth singling out, as it gives rise to the operator 

£ = dd* : R) 4 £ d _ipf ; R) . 

In this case wt = 2 x and Theorem C becomes 

Corollary D (Higher Matrix- Tree Theorem) . For L as above, we have 

det £ = 7x^2 6%. 

T 

Remark 1.6. Variations of Corollary D have appeared in [Ka], [P], 
[DKM] and [L] (note: all but the last reference assume additional con- 
ditions on X, and each work utilizes its own notion of spanning tree). 
When d — 1, we have 9t = 1 = Ox and \ix is the number of vertices of 
X, so we obtain the classical Kirchhoff matrix-tree theorem. Theorem 
C is actually a special case of a more general result, Theorem 6.6 below. 

Corollary D admits the following simpler reformulation (cf. Corol- 
lary 6.10 below). 

Addendum E. 

det£ = ^£ T = ^/i T , 

T T 

where £ T = dd* T : B d ^(T; R) A B d _i(T; R). 

Reidemeister torsion counts spanning trees. Milnor [M] intro- 
duced the notion of Reidemeister torsion t(C*) of a not necessarily 
acyclic finite chain complex C* over a field in which a preferred basis is 
chosen for C* as well as its homology. When C* is the real chain com- 
plex of a finite CW complex X, we will establish a connection between 
the torsion and the enumeration of spanning trees on the skeleta of X. 

Suppose X is a finite, connected CW complex. We give C*(X;R) 
the preferred basis given by its set of cells. We also choose a basis 
for H*(X; R) by selecting a basis for the torsion free part of each inte- 
gral homology group H*(X;1i). Such a basis is called a combinatorial 
basis for the homology and we will denote it by {). The definition of 
Reidemeister torsion r(X; f)) is given in §7 below. 

For k > 0, we define the following quantities: 

• Tk — the set of spanning trees of X^ (for this we require k > 0). 
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• Hk = the square of the covolume of the lattice B k (X;Z) C 
B k (X; R) with respect to the inner product given by restricting 
the standard inner product on C k (X;~R). 

• H k (X; Z)o = the image of the evident homomorphism H k (X; Z) — >• 
H k (X;R). 

• r]k = the square of the covolume of the lattice H k (X;Z) C 
H k (X;M.), where we give Hk(X;M.) the inner product defined 
by identifying the latter with the orthogonal complement of 
B k (X; R) C Z k (X; R) using the inner product arising from the 
standard one on C k (X; R). 

• 9 k = the order of the torsion subgroup of H k (X; Z). 

With respect to the above, we set 



Then 5 k is defined entirely in terms of X. 

Theorem F (Torsion- Tree Theorem). For a finite, connected CW com- 
plex X , we have 

T 2 (x;v=n(5 k y: 4r i)k , 

k>o TeT k+1 

where Q\ denotes the order of of the torsion subgroup of H k (T; Z) for 

re W 

Conventions. We assume the reader is familiar with basic linear algebra 
as well as a first year course on algebraic topology. The topological 
spaces of this paper are equipped with preferred CW structure and 
when we write H*(X; A), we mean cellular homology with coefficients 
in an abelian group A (in practice, A is either Z or 1). If X is a CW 
complex, we write X^ for its fc-skeleton and X k for its set of k cells. 
Thus 

x (k) = x {k-i) u x D hj ^ 

where the union is amalgamated along the attaching map X x S k ~ 1 — > 
X^ -1 '. The k-th Betti number (3 k (X) is defined to be the rank of 
the vector space H k (X;M.). If A is a commutative ring, then the k-th 
real chain group C k (X; A) is by definition the relative homology group 
H k {X^ k \X^ k ~^\ A), which is just the free A-module having basis X k . 

Outline. In §2 we develop basic results about higher dimensional span- 
ning trees. In §3 we prove Theorem A and Addendum B. In §4 we 
prove Theorem C up to identification of the normalizing constant '-fx- 
In §5 we introduce the low temperature limit and use it to show that 
for sufficiently well-behaved W, the determinant of £ tends in the low 
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temperature limit to the determinant of C T , where the latter is defined 
using a spanning tree T in place of X. This result is employed in §6 to 
identify 7x, thereby completing the proof of Theorem C; in so doing 
we generalize Theorem C to Theorem 6.6. Lastly, in §7, we outline 
Milnor's definition of Reidemeister torsion and prove Theorem F. Also, 
in Theorem 7.11, we obtain a different expression for the Reidemeister 
torsion that is expressed in terms of both spanning tree and homology 
truncation data for X. 

Acknowledgements. The authors thank the Los Alamos Center for Non- 
linear Studies and the T-4 division for partially supporting this re- 
search. This material is based upon work supported by the National 
Science Foundation under Grant Nos. CHE-1111350, DMS-0803363 
and DMS-1104355. The third author wishes to thank Paul Kirk for 
discussions related to the beginning of §7 as well as Paul Penfield for 
providing him with a translation to Weyl's article [W]. 

2. Spanning Trees in higher dimensions 

Definition 2.1. Let A be a finite connected CW complex. A /c-cell 
b G Afc is said to be essential if there exists a fc-cycle z G Zk(X;M.) 
such that (z, b) ^ 0. 

Lemma 2.2. Assume in addition X has dimension d. Then adding or 
removing an essential d-cell from X increases or decreases /3d(X) by 
one, respectively, and fixes (3 d _ 1 (X) . 

Proof. Construct a decreasing filtration Y % on X by removing the d- 
cells of X one at a time, X = Y n D Y n ~ x D . . . D Y° = X d ^. Then 
for 1 < j < n, we have an exact sequence in homology 

-> Ha(Y^~ 1 ) H d {Y*) ffd-iQ"'" 1 ) -> H^Yi) -> 

The above factors into two short exact sequences 

-> HaiYi- 1 ) -> H d (Y j ) ^imd^O 

O^Z/imd^ Ha^- 1 ) H d _i(Y j ) ^0, 

where imS* is the image of d*. If the attached cell is essential, then 
im d* is a nontrivial subgroup of Z. Therefore, the first sequence im- 
plies (3 d {Y j ) = P d {Y j ~ l ) + 1, while the second implies p d -i(Y j ) = 
fid-iiY^ 1 ). We may view Y^ as a complex with an additional essen- 
tial cell, or Yi~ l as a complex with an essential cell removed. □ 

Lemma 2.3. X has a spanning tree. 
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Proof. If H d (X; R) = then X is a spanning tree. If H d (X; R) ^ 0, 
then we can pick an essential d-cel\ and remove it, decreasing fid (X) 
by one. Repeat this process until fid is zero. Evidently, the resulting 
subcomplex T contains Xj-i and by Lemma 2.2, we have fid-i(T) = 
j3d-i(X). Hence, T is a spanning tree. □ 

The following is straightforward, and its proof is left to the reader. 

Lemma 2.4. Any spanning tree for X may be obtained by removing 
essential d- cells. Furthermore, ifT is a spanning tree of X, the number 
of essential d-cells withdrawn to construct T is equal to fid(X). 

Lemma 2.5. Let T be a spanning tree of X and let T = Tub, where b 
is an essential cell in T. If b' is an essential d-cell ofT different from 
b, then U := T \ b' is a spanning tree. 

Proof. Since b' is essential, Lemma 2.2 implies Hd(U) has rank zero. 
This lemma also implies fid~i(U) = fid-i(T) = fid-i(T). Since our 
construction leaves the d — 1 skeleton fixed, U is a spanning tree. □ 

Lemma 2.6. Let T be a spanning tree of X and let b G Xd\Td- Then 
[db] generates a torsion element of Hd-\(T] Z). 

Proof. Since T is a spanning tree, we have that b is attached to T 
along its attaching map db — > T. Hence, the homology class [db] lies in 
/f d _i(T;Z). The isomorphism H^T^R) = H d -i(X;R), along with 
the fact that db bounds the cellular chain b in X, implies db is torsion 
inHd-^T-Z). □ 

Recall the linear transformation T: Cd(X;M) —> Zd(X;M.) defined 
in the introduction, which was defined on essential cells as 

f{b) = c/t b , 

where c is a generator of H d (T U 6; Z) and 

h = (c, 6) , 

where the inner product is taken in Cd{X; M.) (here we are using the 
inclusion Hd(TUb; K) C Cd(X; R) to make sense of the inner product). 
Observe that is the order of [db] G Hd-i(T; Z). 

Lemma 2.7. Let T be a spanning tree of X , let bi G Xd \ Td and let 
bj be an essential d-cell. Then 

(T(bi), bj) = ^ . 

Proof, (f (bi), b 3 ) = (c/t bi , b,) = l/t bi (c, b ) = t bj /t bi . □ 
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Corollary 2.8. If U is a spanning tree obtained by adding and then 
removing an essential cell from a spanning tree T as above, then 

(f(b t ),b,)(b l ,U(b J )) = l. 

Lemma 2.9. For an essential d-cell b, the class [db] G H d _i(T;Z) is 
a torsion element of order |t 6 |. In particular, there is a short exact 
sequence 

Z/t b Z -»■ Hd-i(T] Z) H d _±{T U 6; Z) -> . 

Proof. By Lemma 2.6, [96] is a torsion class. Let t be its order. 

By slight abuse of notation, we let db denote the cycle representing 
[db]. Then tdb is also a cycle, which is also the boundary of a unique 
integral d-chain w G C d (T U 6; Z). It is straightforward to check that 
tb — w is a generator of H d (T U 6; Z) = Z^(T U 6; Z). Then (tb — 
w, 6) = t. It follows that t = ±U. The short exact sequence is a direct 
consequence. □ 

For a finite CW complex Y of dimension d, let #y denote the order 
of the torsion subgroup of H d _i(Y; Z). 

Corollary 2.10. For T,U,bi, and bj as above, 

%(T{b i ),b j )=&i J (b i ,U{b j )). 

Proof. Set ti := tj,. and let Y = T U 6, = U U bj. Then the exact 
sequence 

-> Z/tiZ -> F d _i(T; Z) -> H d ^(Y; Z) -> 
gives |ti|0y = ^ and similarly |tj|^y = 6*^/. Consequently, 

e 2 T (f(k), b 3 ) = elutj = efj(b u u(bj)) . □ 

3. Proof of Theorem A and Addendum B 

The proof will proceed along the lines given in [NS] in the classical 
setting. Given a spanning tree T, let {6i, . . . , b k } elements of X d \ T d . 

Lemma 3.1. The collection T(b\), ... ,T(bk) is a basis for Z d (X;M.). 

Proof. Recall that Z d (X; E) = H d (X; R). Let q: X -»■ X/T be the quo- 
tient map. Then the homomorphism : H d (X; M.) —> H d (X/T; R) is an 
isomorphism, and H d (X/T; R) is the vector space with basis 6i, . . . , 6^. 
It's straightforward to check that q*oT: C d (X; R) — > H d (X/T; R) maps 
a d-cell b to itself when 6 G \ T d and is zero otherwise. □ 

Corollary 3.2. For any z G Zd(X;R) ; we /iai>e T(z) = z. 
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Proof. Use the Lemma 3.1 to write z = ^ SiT(bi). Then 

T(z) = SiT(T(bi)) = s ^{h) = z. □ 

i i 

Lemma 3.3. For distinct d- cells hi, bj G X , let %j be the set of span- 
ning trees such that (T(bi),bj) ^ 0. Then 

Y w T {f(bi),bj) R = Y wu{bhU(bj)) R . 
Te%j ueTji 

Proof. From the definition of of the weights, have rjWx = TiWjj. Note 
that (T(bi),bj)n = rj. Using Corollary 2.10, we infer 

6 2 T (T(b l ),b j ) R = 6 2 u (b l ,U(b j )) R . 

Now sum up over all T □ 

Proof of Theorem A. Consider the operator F := Y1t w tT, where the 
sum is over all spanning trees of X. For any pair of rf-cells bj, and bj of 
X we have 

C£2w T f(b i ),b j ) R = ^2 «*r(r(&i),&i)fl 
T Te%j 

= Wu(pi,U(bj)) R by Lemma 3.3 , 

= (b l ,J2^uU(b J )) R 

u 

= {b^^WTfip^R 

T 

Hence F is self-adjoint in the modified inner product. 

If z is a cycle, then F(z) = (52t w t) z =: Az. Consequently, (1/A)F 
restricts to the identity on Zd(X]M). As (1/A)F is self-adjoint, it is 
the orthogonal projection in the modified inner product. □ 

Proof of Addendum B. Let z be the orthogonal projection of i? _1 V in 
the modified inner product. Then i? _1 V - z6 B R (X; M), i.e., 

= (R^V - z, z') R = (V - Rz, z') 

for all z' G Z d (X;R). Hence, V - Rz G B d (X;R). The uniqueness of 
z is a consequence of the fact that B d (X; R) is the orthogonal comple- 
ment to Zd(X; R) in the standard inner product. 
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The proof of the last part is given by direct calculation using the 
self-adjointness of the operator wtT: 

(z,b) = ±(z,b) R , 

T 
T 

4. A WEAK FORM OF THEOREM C 

The goal of this section is to show how Theorem A implies Theorem 
C up to the identification of the prefactor 7. The prefactor will be 
computed in §6, where in addition we prove an enhanced version of 
Theorem C. 

Recall the given function r: X d — > E + of §1. It is convenient to set 

W := lnr: X d — > E . 

Then r b = e w \ and we may also write R = e w : C d (X; E) -»• C d (X; E). 
Conversely, given any function W : X d — > E, we set r := e w : X d — > E + . 
It is convenient to think of W as lying in Cd(X; E) by representing it 

Recall that to each spanning tree T we associated the weight 

where ^ is the order of the torsion subgroup of H d _i(T; Z). 

Remark 4.1. Let M be a smooth manifold and let V be a finite di- 
mensional real vector space. Suppose /: M — > V is a smooth map. 
Then the directional derivative defines a ^-valued, smooth, differential 
1-form df G Q 1 (M;V). In the special case when M — U is a finite 
dimensional real vector space, then Q l (M;V) can be identified with 
the space of smooth maps U — > hom(£7, V). 

Consider the linear operator 

dd* R = de~ w d* : C d ^(X; E) -)■ ^(X; E) . 

Since the image of <9<9|j, is contained in 5^_ 1 (X;E), restriction of this 
operator to B d -i(X;M) gives an isomorphism 

(4) C(W): Bd-i(X;M) A (X; E) . 

For i? = e w , C{W) is the operator C R defined in §1. 
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We can regard W i— > £(W) as defining a smooth map 

(5) £:C d (X;R)->end(B d ^(X;R)), 

which is a family of linear operators parametrized by C d (X;R). To 
avoid notational clutter, when W is understood, we will often write 
£{W) without referring to its argument. Therefore, £ can refer to 
either (5) or (4). 

Proposition 4.2. Theorem A implies the identity 

d In det £ = d In wt ■ 

T 

Remark 4.3. In keeping with our notational ambiguity, the left side of 
the display in Proposition 4.2 is to be interpreted as the value at W of 
d In det £ G Q 1 (C d (X; E); E). 

Proposition 4.2 is equivalent to the statement 

det £ = 7 

T 

for a suitable positive constant 7, as yet to be determined. This gives 
Theorem C modulo the determination of the prefactor 7. 

Proof of Proposition 4-2. We take the differential of the natural loga- 
rithm of det £: 

d In det £ = d tr In £ 

(6) =trd(ln£) 

= tr(£- 1 c/£), 

where d£ = dde~ w d* = -ddWe~ w d*. 
The cyclic property of the trace implies 

(7) tr(£- 1 dC) = -tr(ddWe- w d*C- 1 ). 

If we set A := e~ w d*£^ : B d ^(X; E) E), then tr(£- 1 d£) = 

-tr(<9(iWM) = -tr(dWAd). Consequently, 

dtr\n£ = -ti(dWAd) 

= -J2 (b\dWAd\b) 

bex d 

( 8 ) = - J2 (b\dWA\db) 

b&X d 

= - J2 dW b (b\A\db) , 

bdX d 



KIRCHHOFF'S THEOREMS IN HIGHER DIMENSIONS 13 

where dW b denotes the 6-coordinate function of dW, i.e., dWb(x) = 
dW(x)(b) = W(b), and (i\H\j) stands for the inner product (i,H(j)). 

By definition, A is a left inverse to d: B^(X;R) ->■ B d _i(X; E), so 
the expression (6|^4|9&) is the same as (b,Pb), where P: Cd(X;W) — >■ 
Bji^X] M) is the orthogonal projection in the modified inner product 
( , )r. By Theorem A, we have 

(9) P = /_^^ U ; T f, 

T 

where / is the identity operator. By inserting this expression into 
{b, Pb) and doing some rewriting, we obtain 

(10) (b\A\db) = i £ w T , 

T,b&T d 

where A = wt and the displayed sum is over all spanning trees T 
for which b lies in T. This allows us to rewrite the expression appearing 
in the last line of Eq. (8) as 

(11) £ dw b (b\ A\db) = ^J2Y1 WTdWb ■ 

beX d T b£T d 

On the other hand, for any spanning tree T we have 

(12) d In wt = ~r du>r, 

T T 

where g?u>t is given by 

(13) (Iwt = 4 ^11 e ~ Wb = ~J2 dW b w T . 

b£T d beT d 

Inserting Eq. (13) into Eq. (12) gives 

(14) d In w t = ^ w TdW b . 

T T beT d 

Assembling equations (6), (8), (11), 12 and (14), we conclude 

d\ndetC = -^^2^2w T dW b = dln^2w T . □ 

T beT d T 

5. The low temperature limit 

Here we compute det £ in the low temperature (3 — > oo limit for a 
certain kind of W. We set 

C = de-? w d*: B d ^(X;R) -> B d ^(X;R) 

where W : X d — > M, and j3 G R+ represents inverse temperature. 
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Our freedom in choosing W shows this determinant will tend to 
det £ T , where £ T is £ restricted to a spanning tree. 

Definition 5.1. Fix a spanning tree T of X. A function W: Xd — > R 
is good if 

min W a for any 7 G X d \ T d , 

where k is the number of <i-cells of X. 
Proposition 5.2. For good W : Xd — > ffi, we /icwe 

v det£ T , 
hm — = 1 . 

/?->oo det £ 

Before commencing with the proof, recall the boundary of a d-cell 
a G Xd is given by 

da= ^ b »jJ 

where b a j := (da,j) is the incidence number of a and j. With respect 
to the standard inner product, the adjoint operator d* on a (d— l)-cell 
j is given by 

where b* a := A straightforward computation of the matrix ele- 
ments of £ yields 

£ = V -6 ■ 

a&X d 

Proof of Proposition 5.2. Since Crf(X;M) is a real vector space with 
basis spanned by the set of <i-cells, X^, we have an orthogonal projection 
Q : C d (X; R) -»■ C d (T; M),. This allows us to write 5* = <9£ + 9*, where 
d?p is defined via the commutative diagram 

Ba-i(X; E) >• Cd(X; M) 




It also enables us to write 

£ = £ T + <5£ , 
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where C T = dd^. Together, these imply 

C T = V" e~ ma b b ■ 

a&T d 

Our choice of good W implies any e~^ Wry appearing in the expansion of 
5C must be less than any e~^ Wa appearing in L T and conversely. This 
also means the matrix elements of 5C can be written as a similiar sum; 
the only difference is we instead sum over a G Xd \ 

To simplify taking the limit, we compute the quotient of det C by 
det£ T and let — > oo. Since det C T ^ 0, we may write 

det(£ T + 5C) _ det(J + (C T )- 1 SC) det£ T 
det C T ~ det C T ' 

It suffices to prove that (£ T ) _1 &C tends to the zero operator as (3 — > oo. 
Equivalently, it is enough to show that the matrix elements of (C T )~ 1 5C 
converge to zero. The first bound is of Cjf. 



\C T \ < V e~ pWa \b b 



a€T d 



< e~ PrahlaWa ^\b ai b 



ai'-'aj | • 

This can be further bounded by defining B T = max^ \b a ib a j\. 
Hence, we have 

(15) |4| < e -^in aeTd W a ) B T_ 

The standard formula for the inverse of a matrix gives 

where is the (i,j)-th cofactor of A. Using the exact expression for 
the determinant of C T appearing in Eq. (21) below, Eq. (20) below 3 
and the bound Eq. (15) in the case of the cofactor Cfp we obtain the 
estimate 



«£T% < 



e -^(min aeTd W a ){n-l)/ n _ 



where gx = det(d^dx) depends only on T. 
We bound the elements 5C similarly by 

\8£jk\ < e-P {m[n ^ x ^ T d W "< ) B f 



3 There is no circularity here; Eqs. (20) and (21) do not depend on the material 
in this section. 
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where B T is defined in the obvious fashion. 

Finally, the matrix elements of (C T )~ 1 dC then satisfy the following 
inequality: 

, T ! , (n - l)! e ^(^ Q W Q )(n-l)^T^-l ne -/3min 7 W/ 7jB f 

Collecting terms independent of (3 into AT, we see 

((£ T )~ 1 9£). < e ~' 3 (( n ~ 1 ) min <* W <*-Y, a W Q +min 7 W 7 )yy 

where a 6 and 7 £ A^ \ T^. Our choice of W forces the matrix 
elements to as (3 — > 00. Therefore, 

det£ ^ 
hm - — —7; = det j = 1 , 
/3^oo det£ T 

completing the proof. □ 



6. A GENERALIZED FORM OF THEOREM C 

In this section we will identify the prefactor 7 appearing in Theorem 
6.6. We will also generalize Theorem C in a significant way. 

Covolume. If A is a finitely generated abelian group we let 

A R := A ® z R 

denote its realification, and we let (3(A) = dim^A^ denote the rank of 
A. Let t(A) be the order of the torsion subgroup of A. 

For a homomorphism a : A — > B of abelian groups, we denote 
cck : Am — > £? K be the induced homomorphism of real vector spaces. 

Definition 6.1. A homomorphism a : A — > B of finitely generated 
abelian groups is called a real isomorphism if the induced homomor- 
phism om : — > £?r of real vector spaces is an isomorphism. 

Clearly, a is a real isomorphism if and only if its kernel and its 
cokernel are finite. If a is a real isomorphism, then (3(A) = (3(B), 
where we recall that (3(A) is the rank of A. We will henceforth assume 
that A and B are free abelian. In this case a is a real isomorphism if 
and only if a is a monomorphism with finite cokernel. 

Definition 6.2. For a: A — > B a real isomorphism with A and B free 
abelian, we let 

t(a) £ N 

denote the order of the cokernel, i.e., t(a) := t(B/a(A)). 
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An ordered basis for A determines an ordered basis for A^, and given 
any pair of ordered bases for A, the associated change of basis matrix for 
Ar has determinant ±1. This defines an equivalence relation on ordered 
bases for A with exactly two distinct equivalence classes. A choice of 
equivalence class is referred to as an orientation of A. Consequently, 
when orientations for A and B are chosen, and a: A — > B is a real 
isomorphism, then the determinant det a G R is defined and depends 
only on the choice of orientations. Furthermore, its absolute value 
| det a | is well defined and does not depend on the choice of orientations. 
The latter has the following interpretation: choose an ordered basis for 
B. This defines an inner product on .Br making the ordered basis for 
B into an orthonormal basis for .Br. Then a(A) C -Br is a lattice and 
| det a | is its covolume, that is, the volume of the torus B^/a(A) with 
respect to the induced Riemannian metric, or equivalently, the volume 
of a fundamental domain of the universal covering .Br — > B^/a(A). 

Proposition 6.3. For a real isomorphism a: A — >■ B of finitely gen- 
erated free abelian groups we have | deta| = t(a). 

Proof. Choose an ordered basis for B, and give .Br the induced inner 
product. 

Consider the inclusions a(A) C B C -Br. Then we have a finite 
covering space 

B/a(A) -> B R /a(A) -> B R /B , 

in which the covering projection B R /a(A) — > B^/B is a local isometry 
and B-&/a(A) is the fiber over the basepoint. This shows that the 
covolume of a (A) is the product of the covolume of B C -Br with 
\B/a(A)\ = t(a). But the covolume of B C -Br is 1. □ 

Generalization of Theorem 6.6. Recall that for W : X d — > E, we 

have the operator 

C(W) = de~ w d*: B d -i(X; R) 4 B d ^(X;R) 

which is just C R = dd* R) as defined in the introduction, with R = e w . 
Again, we suppress the argument W from the notation and refer to 
C(W) as C. 

As we showed earlier in Proposition 4.2, we have the following rep- 
resentation: 

(17) det£ = 7^w T , 

T 

where the constant 7 is still to be determined. 
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Definition 6.4. Let A C Cd~i(X; Z) be a subgroup. Define a natural 
number 

H(A) G N 

as follows: let {e^} be a basis for A. Consider the matrix g whose 
(z,j)-entry is given by = (ei,Sj), where the inner product is taken 
in R). Set fi(A) := det g. 

Since e» expressed in the standard basis for Cd_i(X; M.) has integer 
components, we infer that G Z, so is an integer. Alternatively, 
one can define n(A) as the square of the covolume of the lattice A C A®. 
given by restricting the standard inner product of Cd-i(X;M.) to A®.. 
The equivalence of the two definitions can be seen as follows: let B be 
the matrix whose rows are the vectors expressed in an orthonormal 
basis for C^-i- Then | det B\ is the covolume of A C A^. Furthermore, 
g = BB*, so n(A) =detg = (det Bf G N. 

For any abelian group U, we set 

Bd-\ '■= B d ~i(X] U) , 

that is, the image of the boundary operator d: Cd(X; U) — >■ Ca-i(X; U) 
of the cellular chain complex of X with U coefficients. 

Hypothesis 6.5. The inclusion A C Cd-i(X;M.) is such that the or- 
thogonal projection Pa ■ -B*_x ~~ ► A^ is induced by a real isomorphism 
p A : -> A, i.e., P A = (p A ) R . 

Consider the composite operator 

£a ■ A R — > Am 
defined by C A = P A de~ w d*\ All . 

Theorem 6.6 (Generalized Higher Weighted Matrix- Tree Theorem). 
We have 

(18) detC A = lA^2w T , 

T 

where the prefactor is given by 

v(A)t(p A ) 2 



(19) 7 A 



9 



2 



X 

Remark 6.7. The choice A = Bd-i{X; Z) gives Theorem C. 

Remark 6.8. If A = As is the free abelian group generated by a judi- 
ciously chosen subset S C Xj-i, we will obtain fi(A s ) = 1. Using this 
choice of A as well as W = 0, Theorem 6.6 gives a generalization of the 
main result of [P] to CW complexes. 
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Proof of Theorem 6. 6. As above, we have 

C := dd* R = de~ w d* : B d ^(X; R) A B d ^(X; R) . 

Then 

I~-a = Pa£Pa j 

which implies 

(20) det C A = det(£) Aei(P A P* A ) . 

If we apply this to Eq. (17), we reproduce Eq. (18) with 7a = 7 det (Pa-Pa)- 
It suffices to identify the prefactor 7a- 

Consider the operator C T for some spanning tree. We have 

(21) det C T = det{d T e~ w d* T ) = det{d T d* T e~ w ) = f£ det(<9 T <9* ) 

Applying Eq. (18) in the case of good W and in the low temperature 
limit, the left hand side of that equation tends to the determinant of the 
operator C7 A for the spanning tree Tclof maximal weight, whereas 
the right hand side is dominated by the single contribution associated 
with the same spanning tree T. Consequently, Eq. (21) implies 

(22) det (d T d* T ) det (Pj(Pl)*) = 1 A 6 2 T . 

Since Pj = (pa)»> where the real isomorphism p T A : B d _i(T; Z)) — > A 
is obtained by composing the real isomorphism p A '■ Pd_i(X;Z) — > A 
with the inclusion B d _i(T; Z) C B d _i(X; Z), we have 

det (P A P* A ) = fi(A)(fi(B d ^(T; Z)))~\detp T A ) 2 . 

We further note that, since T is a spanning tree, the free abelian group 
B d _i(T; Z) has basis {ctei, . . . , de s }, where e±, . . . , e s are the fc-cells of 
T, so that we have a matrix g of inner products with the matrix ele- 
ments gij = (drei, drej) = (d^dxei, e,), which implies /i(Pd_i(T; Z)) = 
det(d^dT) ■ Then Eq. (22) assumes the form 

fx(A)(detp T A ) 2 = lA 9 2 T . 

Combining this with Proposition 6.3 results in 

l 2S ) 1A = ^ • 

Urp 

The right side of Eq. (23) is written in terms of a particular spanning 
tree T, however, it does not actually depend on this choice. An in- 
variant expression that does not contain T is obtained by using the 
following relations: 



(24) ^4 = ^iC*; niB^{P; Z)) = 

t{pA> X 
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Substituting Eq. (24) into Eq. (23) results in an invariant expression 
for 7^4, given by Eq. (19). □ 

Alternative forms of Theorem C. In this subsection we deduce 
Addendum E as well as a generalization of it to the weighted case. Let 
us now return to the more general situation of Theorem 6.6. 

Theorem 6.9. With A C Cd-i(X;Z) as above, we have 

det La = det L T A . 

T 

Proof. Using Eq. (24) we infer that 

^A)t{P A f ^A)t(P T A f 

for any spanning tree T. Combining this with Theorem 6.6 in the case 
of a spanning tree T we obtain 

det C T A = 'JaWt . 

The conclusion now follows by summing over all T. □ 
In the special case when A = B^i{X; Z), Theorem 6.9 reduces to 
Corollary 6.10. det £ = ^T^ et ^ T = ^t^t- 

7. Reidemeister torsion and Theorem F 

Reidemeister torsion. Milnor [M] defined the Reidemeister torsion 
of a not necessarily acyclic finite chain complex over a field equipped 
with the auxiliary structure of an ordered basis of its chains as well 
as a choice of ordered basis of its homology groups. In this section we 
restrict ourselves to the case of torsion for chain complexes defined over 
the real numbers. 

Consider the case of a chain complex C* of finite dimensional vector 
spaces over K having non-trivial terms in degrees < * < d. Let 
d: Ck — > Ck-i be the boundary operator. Let C Ck be the subspace 
of /c-cycles and let P>k C Zk the subspace of /c-boundaries. We also set 
Hk = Z k /B k . 

We then have short exact sequences 

-> Z k -> C k -> B k -i ^0 and ->■ B k -> Z k -> H k -> . 

If we choose splittings s^-i '■ B^-i — > Ck and t^' Hk — > Zk, we are 
entitled to write C k = Z k ® B k = B k ®H k ®B k . 
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Pick bases b k := {b\},c k := {c l J, f) ft := {h\} for B k ,C k , and if fc , 
respectively. It follows that {sfc_i(6^._ 1 ), t k (h k ), b l k } forms another basis 
for C k . Let {b k \) k b k _i\ denote this basis and let 

[Mfcbfc-iAk] 

denote the change of basis matrix that expresses the basis b k h, k b k -i in 
terms of the basis c^. Let c = {c k } and f) = {b, k }- 

Definition 7.1 (Milnor [M, p. 365]). The torsion of the pair (C*, f)) is 
defined by 

r(a) = n det M*Wcfc] ( " 1)h , 

fc>0 

which is consistent with Milnor's definition with respect to the identi- 
fication of A'i(R) = R x given by the determinant function. 

Milnor shows that the definition is independent of the choice of b 
as well as the splittings. Thus the torsion is really an invariant of the 
triple (C*, c, f)). 

In what follows, = C*(X;R) is the cellular chain complex of a 
finite connected CW complex X which has a preferred basis consisting 
of the set of cells. In this case, we think of the torsion as an invariant 
of the pair (X, [)) and we set 

r(X;b) :=t(C*(X;R)), 

where we have indicated in the notation the dependence on the choice 
of homology basis. It will be useful to single out a specific kind of 
homology basis. Let if„,(A;Z)o C H*(X,M) be the lattice given by 
taking the image of the evident homomorphism H*(X] Z) — > H*(X; M.). 
Note that if*(X;Z) has a preferred isomorphism to the torsion free 
part of Hz(X; Z). 

Definition 7.2. A combinatorial basis for H*(X\ M.) consists of a basis 
for H k (X;Z) for k > 0. 

Henceforth we fix a combinatorial basis h, . Let r : JI k X k — > M + 
be a positive- valued function on the set of cells of X. As in previous 
sections, we write R: C*(X;R) — > C*(X;R) for the linear transforma- 
tion determined by 6 i — r^b, and R = e w . We have a modified inner 
product (b,b')R = (rbb,b'). We also have an operator 

C k (W) = dd* R := de- Wk ^d*e Wk : B k (X; R) ^ B k (X; R) , 

where d* R is the formal adjoint to d: Cfc + i(A;R) — > Cfc(A;R) in the 
modified inner product on both source and target. We define H^X; R) 
to be the orthogonal compliment of B k (X; R) in Z k (X; R) with respect 
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to modified inner product on C k (X; R), and we then have a preferred 
identification Hj}(X;~R) = Hk(X;M) given by sending a cycle to its 
homology class. As in the introduction, we let r] k be the square of the 
covolume of H k (X;Z*) C H k (X;M), with respect to the basis t) k for 
H k (X; Z) and the inner product on Hj?(X; R) obtained by restricting 
the modified inner product on C k (X; IR). 4 

Theorem 7.3. Let X be a finite, connected CW complex. Then 

T 2 (X h) = even det £fc(WQ Qfc odd,b&X k EL even Vk 

cvcn.feGXfe e kb 1 Ifc odd Vk 

Remark 7.4. If we take W = 0, then Theorem 7.3 immediately implies 
that t 2 (X; fj) is an invariant of the lattice H*(X; Z) C H*(X\ R) rather 
than just an invariant of the specific choice of combinatorial basis [). 
Since this lattice doesn't depend on any choices, we infer that t 2 (X; f)) 
depends only on the CW structure of X. In fact, the method of proof 
of [M, th. 7.2] shows that t 2 (X; f)) is invariant under subdivision. 

Proof of Theorem 7.3. For the purpose of this proof we suppress W 
and write £ = C{W). We also set C* := C*(X; R). Define the splitting 
maps s fc _i : B k _ x ->■ C fc by 

Let R) denote the image of Sk-i, and similarly we define B R (X; Z) 

to be Sk-i(B k (X;Z)). Note that B R (X;M) is the orthogonal compli- 
ment to Z k in the modified inner product on C k . 

Let 7 fc denote the square of the covolume of B R (X; Z) C B R (X; IR), 
using the inner product on B R (X; R) induced by the modified inner 
product on C k . Similarly, let 7*._i denote the square of the covol- 
ume of Bk-i(X; Z) C Bk-i(X; R), where B k -i(X; R) is given the inner 
product by restricting the modified inner product on Ch-i- Using the 
isomorphism B k © H k © B k _i -H- C k determined by the splitting, we 
infer 

k 

(25) det[b fc h fc b fc _!/c fc ] 2 - 



so the square of the Reidemeister torsion is 

Ilk even 1kVkl k Ilk odd,beX k 



(26) r\X- h) 



even,b£X k 



4 In the introduction, rj k was denned only in the case when W = 0; the current 
notation applies to an arbitrary W. 
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Since Sk = d^rC k , we have s* k = C k l d (since C is self-adjoint). There- 
fore, s* k Sk = C^dd^C^ 1 = C^ 1 . We use this fact to compute the 
quotient of 7 fc /ju-x- Recall that 7 fc is given as the determinant of the 
inner product matrix, we compute 

(s fc - 1 (6fc_i)»Sfc-i(6fc- 1 ))fl = (Sfc-i5fc-i(6fc-i),(6 , fc-i))fl 

= {ct-M-i)M-x)) R - 

The determinant of the matrix with these entries latter is, by definition, 
(det U) 2 det where U is the change of basis matrix expressing 

bk-i in terms of an orthornormal basis for Bk-i(X; R) in the modified 
inner product. A similar observation shows that the determinant of 
the matrix whose entries are (b l k _ 1 ,b k _ 1 )R is (detf/) 2 , and this is just 
7fc-i- 

Consequently, the quotient of these determinants is 

7 fc 1 

(27) 



7 fc _i det £ fe _i ' 

Inserting Eqn. (27) into Eqn. (26) and performing the evident cancel- 
lations, we conclude 

T 2 (X- h) - even^ et ^k(W) fife odd,fegX fc EL even Vk m 

V IL odd det C k (W) Uk even,beX k e ** 1 Lk odd Vk 

In the special case when W = 0, we can combine Theorem 7.3 with 
Corollary D. This immediately gives Theorem F: 

Corollary 7.5 (Torsion- Tree Theorem). For a finite, connected CW 
complex X , we have 



r 2 (X;h) = n(4 a 

k>0 TeT k 



2\(-l)» 

t) 



An alternative formula. In this part we shall derive a different for- 
mula for the torsion in terms of a single spanning tree in each degree 
as well as a choice of auxiliary structure, namely, homology truncation 
data for X. 

Hypothesis 7.6. For each k > 1, we fix a spanning tree T k for X™ . 
Our convention is to set T° = 0. 

Definition 7.7. A homology truncation of X in degree k > 0, subor- 
dinate to T k , is a subcomplex i: V k C X^ such that T fc C V k and 
i*: H*(y ;R) — >■ if*(AT;IR) are isomorphisms for * < k. 
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In induction argument similar to the proof of Lemma 2.3 shows that 
homology truncations exist. Note that V° consists of a single vertex of 
X. We have a filtration 

T° C V° C • • • C X {k ~ 1] C T k C^C X (k) C • • • 

Lemma 7.8. The choice of spanning tree T k determines a splitting 
Bk-i(X;W.) — y C k (X;R). The choice of homology truncation V k sub- 
ordinate to T k determines a splitting H k (X;R) —y Z k (X;R). 

Proof. The first splitting is the composition 

Bk-i(X; R) = B k _i(T h ; R) B k (T k ; R) -y C k (X; R) 
and the second is given by 

H k (X; R) ^y H k (V k ; R) = Z k (V k ; R) ^ Z k (X; R) . □ 

Define a basis for _B fc (X;Z), b fc = {6*}, as given by the cells of T k . 
Here we are using the preferred isomorphism B k (X;Ij) = C k {T k ;'L). 
This defines a basis b k -i for B k _i(X; R) by = 96^}. The basis for 
homology in degree is the combinatorial basis i) k given as an input 
to the torsion. As always, the basis for C k (X; R) is given by the set of 
fc-cells. 

Before explicitly identifying the torsion, note that in each dimension 
k there are essentially three types of cells: 

X k = (T k ) U (V k \ T k ) U (X k \ V k ) . 

Roughly speaking, the first set of cells contributes to B k , the second 
set contributes to H k and the last set contributes to B k . This gives us 
a decomposition of the fc-chains 

(28) C k (X; R) = C k (T k ; R) © C k (V k /T k ; R) © C k (X/V k ; R) . 

(when k = 0, we replace C (V°/T°;R) with C (V ,T°;R) = R, etc.) 

We first identify the homological contribution in degree k to the 
torsion. With respect to the splitting Eq. (28), the combinatorial basis 
f)fc has image contained in the direct sum 

C k (T k ; R) © C k (V k /T k ; R) = C k (V k ; R) . 

Hence, its contribution to the torsion is left invariant if we project 
these elements onto C k (V k /T k ;R) = H k (V k /T k ;R) = H k (X;R) (since 
the other summand C k {T k ;R) = B k (T k ; R) maps to B k (X;R) and 
the relevant determinant remains unchanged if we project away from 
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B k (X;R)). Consequently, the homological contribution to the torsion 
in degree k is given by the determinant of the composite 

H k (X; R) ^ H k (V k - R) ^ H k (V k /T k ; R) , 

where p: V k — > V k /T k is the quotient map. So we wish to identify 
det p*/ det z*. 

Definition 7.9. Let 

XkEN 

denote the square of the determinant of z*: H k (V k ;R) — > H k (X;R), 
i.e., the square of the covolume of the lattice i*(H k (V k ; Z)) C H k (X;R). 

Applying Proposition 6.3 to the real isomorphism H k {V k ] r L) — > 
H k (V k /T k ;Z), we infer 

Lemma 7.10. The determinant of p* is the ratio ±6 T k /9 V k . 

Consequently, up to sign, the contribution of fa to the determinant 
defining the Reidemeister torsion is 

(29) 

Vv k V Xk 

We next identify the contribution in degree k to the torsion provided 
by the basis b&. As defined above this basis is given by the boundaries 
of the cells of Tfc+i- This leads us to consider the composite 

(30) C k+1 (T k+1 ; Z) A B k (T k+1 ; Z) A C k (X/V k ; Z) , 

where q k is induced by the quotient map T k+1 — > X/V k . The ho- 
mo morphism d is an isomorphism and so it has determinant ±1. The 
second homomorphism q k is a real isomorphism and therefore the deter- 
minant of its realification, det((gfc)iR), has value ±t(q k ) by Proposition 
6.3. Note that (qk)w. is the restriction of the orthogonal projection 
C k (X;R) C k {X/V k ;R) to the subspace B k (T k+1 ;R) C C k (X;R). 
and the projection of b k onto this summand gives its contribution to 
the torsion. Hence, the determinant of the composition (q k )u ° d is 
±t(q k ). So the contribution in degree k of b k to the torsion is ±t(q k ). 

Lastly, the contribution to the torsion in degree k provided by the 
basis bfc_i is given by the standard basis of C k (T k ; R) via the splitting 
Eq. (28). It is then evident that the contribution in degree k of b k _i 
to the torsion is 1. 

Assembling, we obtain 

(31) det[b fc f) fc b fc _i] = ±t(q k ) ■ — • 1 . 

Vv k \/Xk 
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Forming the square of the Reidemeister torsion, we conclude 

Theorem 7.11. For a connected, finite CW complex X with com- 
binatorial homology basis f), spanning tree data {T k } and homology 
truncation data {V h }, we have 

where q k : B k (T k+1 ; Z) -> C k (X/V k ; Z) and »eN are as above. 

Example 7.12. If X has dimension one, then all terms appearing in 
Theorem 7.11 are equal to one. Hence, t 2 (X; h) = 1 whenever X is a 
connected finite graph. 

Example 7.13. Let X = M.P 2 . Then we may choose T 1 = * = V 1 and 
T 2 = M.P 2 = V 2 . In this instance, the only non-trivial term appearing 
in Theorem 7.11 is t(q x ) 2 = 4. Hence t 2 (RP 2 ; h) = \. 
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